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Abstract 



We derive explicit distance bounds for Stratonovich iterated integrals along two Gaus- 
sian processes (also known as signatures of Gaussian rough paths) based on the regularity 
assumption of their covariance functions. Similar estimates have been obtained recently in 
[Friz- Riedel, AIHP, to appear]. One advantage of our argument is that we obtain the bound 
for the third level iterated integrals merely based on the first two levels, and this reflects the 
intrinsic nature of rough paths. Our estimates are sharp when both covariance functions 
^vq ' have finite 1-variation, which includes a large class of Gaussian processes. 

Two applications of our estimates are discussed. The first one gives the a.s. conver- 
gence rates for approximated solutions to rough differential equations driven by Gaussian 
processes. In the second example, we show how to recover the optimal time regularity for 
solutions of some rough SPDEs. 
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The intersection between rough path theory and Gaussian processes has been an active 
research area in recent years ( [FVlOaj . |FV10bj . [Hllj ). The central idea of rough paths, as 
realized by Lyons ( |L98j ) . is that the key properties needed for defining integration against 
an irregular path do not only come from the path itself, but from the path together with a 
sequence of iterated integrals along the path, namely 

X? >t = f dX Ul ®...®dX Un . (1) 



X 
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In particular, Lyons extension theorem shows that for paths of finite p-variation, the first 
[p\ levels iterated integrals determine all higher levels. For instance, if p = 1, the path 
has bounded variation and the higher iterated integrals coincide with the usual Riemann- 
Stieltjes integrals. However, for p > 2, this is not true anymore and one has to say what the 
second (and possibly higher) order iterated integrals should be before they determine the 
whole rough path. 

Lyons and Zeitouni ( |LZ98j ) were the first to study iterated Wiener integrals in the sense 
of rough paths. They provide sharp exponential bounds on the iterated integrals of all levels 
by controlling the variation norm of the Levy area. The case of more general Gaussian 
processes were studied by Friz and Victoir in |FV10aj and [FVlObj . They showed that if 
A" is a Gaussian process with covariance of finite p- variation for some p £ [1,2), then its 
iterated integrals in the sense of ([1} can be defined in a natural way and we can lift X to a 
Gaussian rough path X. 

In the recent work |FR12j , Friz and the first author compared the two lift maps X and Y 
for the joint process {X, Y). It was shown that their average distance in rough paths topology 
can be controlled by the value sup t \X t — Y t \^ L2 for some £ > 0, and a sharp quantitative 
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estimate for ( was given. In particular, it was shown that considering both rough paths in 
a larger rough paths space (and therefore in a different topology) allows for larger choices 
of (. Using this, the authors derived essentially optimal convergence rates for X £ — > X in 
rough paths topology when e — > where X c is a suitable approximation of X. 

In order to prove this result, sharp estimates of |X™ t — Y™ t \ need to be calculated on 
every level n. Under the assumption p £ [1, |), the sample paths of X and Y are p- rough 
paths for any p > 2p, hence we can always choose p < 3 and therefore the first two levels 
determine the entire rough path. Lyons' continuity theorem then suggests that one only 
needs to give sharp estimates on level 1 and 2; the estimates on the higher levels can be 
obtained from the lower levels through induction. On the other hand, interestingly, one 
additional level was estimated "by hand" in }FR121 before performing the induction. To 
understand the necessity of computing this additional term, let us note from |L98) that 
the standard distance for two deterministic p-rough paths takes the form of the smallest 
constant C n such that 

|Xj it -Y? it |<C„ea;(M)S, n = l,---,\jp\ 

holds for all s < t where w is a control function to be defined later. The exponent on the 
control for the next level is expected to be 

n + l_bJ+l >l5 (2) 



p p 

so when one repeats Young's trick of dropping points in the induction argument (the key 
idea of the extension theorem), condition @ will ensure that one can establish a maximal 
inequality for the next level. However, in the current problem where Gaussian randomness 
is involved, the L 2 distance for the first [2p\ iterated integrals takes the form 

IX^-Y™^ <C n eLo{s,t)M^, „=1 )2 , pe[l,|), 

where 7 might be much larger than p. Thus, the 'n — 1' in the exponent leaves condition 
^ unsatisfied, and one needs to compute the third level by hand before starting induction 
on n. 

In this article, we resolve the difficulty by moving part of e to fill in the gap in the control 
so that the exponent for the third level control reaches 1. In this way, we obtain the third 
level estimate merely based on the first two levels, and it takes the form 

\Xl t -Yl t \ L i <C 3 eMM), 

where 77 £ (Ojl]> an d its exact value depends on 7 and p. We see that there is a 1 — i] 
reduction in the exponent of e, which is due to the fact that it is used to compensate the 
control exponent. This interplay between the 'rate' exponent and the control exponent can 
be viewed as an analogy to the relationship between time and space regularities for solutions 
to SPDEs. We will make the above heuristic argument rigorous in section 4. We also refer 
to the recent work [LXllj for the situation of deterministic rough paths. 
Our main theorem is the following. 



Theorem 1.1. Let (X, Y) = (X 1 , Y 1 , • • • ,X d ,Y d ) : [0,T] -> R d+d be a centered Gaussian 
process on a probability space (Q,IF,P) where [X l ,Y l ) and (X 3 ,Y J ) are independent for 
i 7^ j, with continuous sample paths and covariance function R(x,Y)'- [0, T] 2 1— > Jj2dx2d^ 
Assume further that there is a p £ [1, |) such that the p-variation of R(x,Y) * s bounded by a 
finite constant K . Let 7 > p such that - + - > 1. Then, for every a > 2j, N > [a\. q > 1 
and every S > small enough, there exists a constant C = C(p,"/,o~,K,5,q) such that 

ft) I f^ + 1 P > 1. then 

\Q»_ var (X,Y)\ L « <C sup iXt-Ytf'*. 

te[o,T] 
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(ii) //i + I < 1, then 

\g?_ var (X,Y)\ Lq <C sup \X t -Y t fe 
te[o,T] 



The proof of this theorem will be postponed to section 4.2 after we have established all 
the estimates needed. We first give two remarks. 

Remark 1.2. We emphasize that the constant C in the above theorem depends on the 
process (X, Y) only through the parameters p and K. 



Remark 1.3. For N — [a\, Q„_ var denotes an inhomogeneous rough paths metric. See 
section^ for the precise definition. 

Our paper is structured as follows. In section 2, we provide some important concepts 
and notations from rough path theory that are necessary for our problem. In section 3, 
we introduce the class of Gaussian processes which possess a lift to Gaussian rough paths 
and estimate the difference of two Gaussian rough paths on level one and two. Section 
4 is devoted to the proof of the main theorem. We first obtain the third level estimate 
directly from the first two levels, which requires a technical extension of Lyons' continuity 
theorem, and justify the heuristic argument above rigorously. All higher level estimates are 
then obtained with the induction procedure in |L98] . and the claim of the main theorem 
follows. In section 5, we give two applications of our main theorem. The first one deals 
with convergence rates for Wong-Zakai approximations in the context of rough differential 
equations. The second example shows how to derive optimal time regularity for the solution 
of a modified stochastic heat equation seen as an evolution in rough paths space. 

Notations. Throughout the paper, C,C n ,C n (p,j) will denote constants depending on 
certain parameters only, and their actual values may change from line to line. 

Acknowledgements. We wish to thank our advisors, Peter Friz and Terry Lyons, 
for their helpful discussions and support during the project. S.Riedel is supported by a 
PhD scholarship from the Berlin Mathematical School (BMS). W.Xu is supported by the 
Oxford-Man Institute Scholarship. 



2 Elements from Rough path theory 

In this section, we introduce the concepts and definitions from rough path theory that are 
necessary for our current application. For a detailed account of the theory, we refer readers 
to |FV10bj . |LCL06] and \LQ02\ . 

Fix the time interval [0, T]. For all s < t £ [0,T], let A s t denote the simplex 

{(ui,u 2 ) | s <Ul <u 2 <t}, 

and we simply write A for Ao,t- In what follows, we will use x to denote an Revalued 
path, and X to denote a stochastic process in R d , which is a Gaussian process in the current 
paper. For any integer N, let 

T N (R d ) = R M. d © ■ • • © (R d )® N 

denote the truncated tensor algebra. The space of all continuous bounded variation paths 
x: [0, T] -> R d is denoted by C 1 ~ var (R d ). For a path x e C 1 ~ var {R d ) 1 we use the bold letter 
X to denote its n-th level iterated tensor integral: 

X" t = / dx Ul (g) • • • <gi dx Un . 

J s<ui<---<u n <t 
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The lift map Sn taking i to a T Ar (R d )-valucd path is dehned by 

N 
n=l 

For a path x, write x Sj t — x t — x s and we have x s j = X s t . It is well known that Sn(x) is 
a multiplicative functional, that is, for any s < u < t, we have 

Sn(x)s,u ® S N (x) u j = S N (x) s . t , 

where the multiplication takes place in T N (R d ). 

For each subspace (R d )®", there is an associated tensor norm | • |. If X,Y are two 
multiplicative functionals in T N (R d ), then for each p > 1, we define their p- variation distance 

by 

e£ oaP (X,Y) :=maocsup(^|X^ ti+l -Y» ti+ J^)t, 

i 

where the supremum is taken over all finite partitions of the interval [0, T] . If N = [p\ , this 
defines 3, rough, paths metric cLnd we only write Qp—var 

(X, Y) in this case. 

Remark 2.1. Note that we only consider so-called inhomogeneous rough paths metrics in 
this paper. The reason for this is that the Ito-Lyons solution map for rough paths is locally 
Lipschitz with respect to these metrics (cf. \FV10b[ Chapter 10]). 

We define the subset G N (R d ) C T N (R d ) to be 

G N (R d ) = {S N (x) ,i : x e C^ var (R d )}. 

The multiplicativity of Sn implies that G N (R d ) is a group with multiplication ® and identity 
element 1. If a; € C 1 - 1 " 1 ''^), one can shows that actually S N (x) s ^ t G G JV (R d ) for all s < t. 

Definition 2.2. A function uj : A — > R + is called a control if it is continuous, vanishes on 
the diagonal, and is superadditive in the sense that for any s < u < t, we have 

u(s, u) + oj(u, t) < oj(s, t). 

We say a multiplicative functional X in T N (R d ) has finite p- variation (p > 1) controlled 
by uj if for each n < N, there exists a constant C n such that for all s < t, we have 

\X? tt \<C n uj(s,t)$. 

Definition 2.3. Let p > 1. A geometric p-rough path is a continuous path in GL p J 
which is in the p-vavicitioTi closure (w.v.t the. ine-tric Qp—var 

) of the set of bounded variation 
paths. We use C°' p ~ var ([0, T], G^J (R d )) to denote the space of geometric p-rough paths. 

By Lyon's extension theorem (cf. [L981 Theorem 2.2.1]), every (geometric) p-rough path 
X can be lifted to a g-rough path for every q > p. Abusing notation, we will use the same 
letter X to denote this lift. We will write X" = 7T„(X) where 7r„ denotes the projection of 
T N (M. d ) onto the n-th tensor product, n < N. If x = 7Ti(X), we will also use the notation 

X" t = / dx® ■■ -®dx 
(even though this integral does not have to exist as a limit of Riemann sums). 



4 



3 2D variation and Gaussian rough paths 



If I = [a, b] is an interval, a dissection of I is a finite subset of points of the form {a = tg < 
. . . < t m — b}. The family of all dissections of I is denoted by T>(I). 

Let I C R be an interval and A = [a, b] x [c, d] C I x I be a rectangle. If /: I x 7 — >■ V 
is a function, mapping into a normed vector space V , we define the rectangular increment 
f(A) by setting 

/<">=-/( ji)=-/(5)-/(;)-/(S) + /(;). 



Definition 3.1. Let p > 1 and f:IxI—>V. For [s,t] x [u,u] C I X I, set 



V P (f; [s,t] x [n,u]) := 



sup ^ 

(i 4 )ex»([a,t]) t t / 
\(t$)ex>(Ku]) " ' 



J 



IfV p (f,I x I) < oo, we saj/ t/iai / /ias finite (2D) p-variation. We also defir, 



Voo(f; [s,t] x [u,v]) := sup 

<T,re[s,t] 



<T,T 



Lemma 3.2. Let f : I x I —> V be a continuous map and 1 < p < p' < oo. Assume that f 
has finite p-variation. Then for every [s, t] x [u. v] C I x I we have 



x Ml) < Voo(f;[s,t] x [u^]) 1 $V p (f;[s,t] x [u,v])* 
Proof. Let (U) e £>([s,i]) and (t'j) G Then, 



E 



<V«, (/;[*,<] x 



Taking the supremum over all partitions gives the claim. 



/ 



ti, ti+l 



□ 



Lemma 3.3. Let f : I x L — > R &e continuous with finite p-variation. Choose p' such that 
p' > P */P — 1 arl ^ p' > p if p > 1. Then there is a control U) and a constant C = C(p,p / ) 
such that 

V pl (f;JxJ)<u(J)V <CV p (f;JxJ) 
holds for every interval J C I . 

Proof. Follows from [FV1JJ Theorem 1]. □ 



Let X = (X 1 , X d ) : L -> R d be a centered, stochastic process. Then the covariance 
function Rx{s,t) := Covx(s,t) — E(X S ® X t ) is a map i?x : I x I ~R dxd and we can 
ask for its p-variation (we will use the letter p instead of p in this context). Clearly, Rx 
has finite p-variation if and only if for every i, j e {1, . . . , d} the map s, t h-» E(X l s Xf) has 
finite p- variation. In particular, if X 1 and X^ are independent for i ^ j, Rx has finite 
p-variation if and only if has finite p-variation for every i = l,...,d. In the next 
example, we calculate the p-variation for the covariances of some well-known real valued 
Gaussian processes. In particular, we will see that many interesting Gaussian processes 
have a covariance of finite 1-variation. 
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Example 3.4 

(i) Let X — B be a Brownian motion. Then Rs(s,t) = min{s,i} and thus, for 

A = [s, t] x [u, v], 

\R(A)\ = \{s,t)D(u,v)\= 5 x=y dxdy. 

J [s,t] x [u,v] 

This shows that Rb has finite 1-variation on any interval /. 

(ii) More generally, let /: [0,T] —> R be a left-continuous, locally bounded function. Set 

ft 



X t = f f(r)dB r . 
Jo 

Then, for A = [s, t] (1 [it, v] we have by the Ito isometry, 



Rx{A) = E 



fdB / fdB 

s,t] J[u,v] 



&x=yf(x)f(v) dx dy 

s,i] X [u,v] 



which shows that Rx has finite 1-variation. 

(iii) Let X be an Ornstein-Uhlenbeck process, i.e. X is the solution of the SDE 

dX t = -8X t dt + adB t (3) 

for some 9, a > 0. If we claim that Xq = 0, one can show that X is centered, Gaussian 
and a direct calculation shows that the covariance of X has finite 1-variation on any 
interval [0,T]. The same is true considering the stationary solution of (J3J instead. 

(iv) If X is a continuous Gaussian martingale, it can be written as a time-changed 
Brownian motion. Since the p-variation of its covariance is invariant under 
time-change, X has again a covariance of finite 1-variation. 

(v) If X : [0, T] — > R is centered Gaussian with Xq = 0, we can define a Gaussian bridge 

by 

-^Bridge (t) = X t — t~^T- 

One can easily show that if the covariance of X has finite p-variation, the same is 
true for ^Bridge- In particular, Brownian bridges have finite 1-variation. 



Next, we cite the fundamental existence result about Gaussian rough paths. For a proof, 
cf. [FVlOaj or [FVlObl Chapter 15]. 

Theorem 3.5 (Friz, Victoir). Let X : [0,T] — > M. d be a centered Gaussian process with 
continuous sample paths and independent components. Assume that there is a p G [1,2) 
such that V p (Rx] [0,T] 2 ) < oo. Then X admits a lift X to a process whose sample paths are 
geometric p-rough paths for any p > 2p, i.e. with sample paths in C°' p ~ var ({0, T], CL?>J (R d )) 
and 7Ti(X Si t) — X t — X s for any s < t. 

In the next proposition, we give an upper L 2 -estimate for the difference of two Gaussian 
rough paths on the first two levels. 

Proposition 3.6. Let (X, Y) = (X^Y 1 ,--- ,X d ,Y d ) : [0,T] -> R d+d be a centered Gaus- 
sian process with continuous sample paths where (X l ,Y l ) and (X J ,Y 3 ) are independent for 
i ^ j. Let p G [1, 1) and assume that V p '{R(x,Y)-> [0,T] 2 ) < K < +oo for a constant K > 
where p' < p in the case p > 1 and p' = 1 in the case p — 1 . Let 7 > p such that ~ + ^ > 1 . 
Then there are constants Cq,C\,C2 dependending on p^p',^/ and LC and a control u such 
that u)(0,T) < C and 

\X a , t - Y s , t \ L2 < Ci sup \X U -Yul^i uj(s,t)^ 
ue\s,t] 
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and 



X SyU ®dX u 



Y s ,u ® dY u 



<C 2 sup \X U - Y U \ L2 1 u(s,t)*y + 2" 
l 2 ue[s,t] 



hold for every s < t. 

Proof. Note first that, by assumption on V P '{R(x,y)\ [0, T] 2 ) , Lemma [3.31 guarantees that 
there is a control to and a constant c\ — c\(p, p') such that 

V P (R X ; [s,t] 2 ) V V P {R Y ; [s,t] 2 )\JV p (R {x _ Y) - [s,t] 2 ) < u{s,t) x 'P 

holds for all s < t and i = 1, . . . , d with the property that u>(0, T) < c\K p —: Cq. We will 
estimate both levels componentwise. We start with the first level. Let i £ {1, . . . , d}. Then, 



XI 



S,t\ L 2 



R. 



S, t 

s, t 



< 



V^(R {x _ Y y,[s,t] 2 ) 



and thus 



\X s , t - Y s>t \ L2 < c 2 ^V^R ix ^ Y) ; [s,t] 2 ). 

For the second level, consider first the case i = j. We have, using that (X,Y) is Gaussian 
and that we are dealing with geometric rough paths, 



XI u dXI 



Y l dY % 

1 s.u u 1 u 



L 2 



\\{Xl t ) 2 ~{Yl t f\ 



L 2 



L 2 



<c 3 \Xi t -Yl t \ T A\XiA L * + \Yl t \ L z) 



From the first part, wc know that 



\K.t - Y s \ t \ L2 < ^JVj(R(x-y)'i [ s ,t] 2 ). 



Furthermore, 



\xIAl 2 



Rx 



s, t 
s, t 



< Jv p {Rx-AsA 2 )<u(s,t)-> 



and the same holds for |Y s l t |i2. Hence 



For i ^j, 



XI „, dX 



s.u u 



Yi u dY: 



< 



L 2 



c 4 ^V 7 (R ( x-y); [s,t] 2 )uj(s,t) 2 p. 



XldXl- Yl „, dYJ 



s.u u 



s.u u 



< 



(X i -Y% tU dXi 



L 2 



Yl u d(Xi-Yi\ 



L 2 



We estimate the first term. From independence, 
E 



t * 2 

{X 1 -Y l ) s .udXi 



s,t] 2 



R(X i -Y i ) 



s, u 

S, V 



dR X i (u, v) 
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where the integral on the right is a 2D Young integral!]] By a 2D Young estimate (cf. |T02j ). 



R(X'-Yi) [ I'y ) dR X i(u,v) 



< c 5 (p,7)^ 7 (^(xi-ri); [s,t] 2 )V p (R xi ;[s,t] 2 ) 

<c a V 1 {R (X -Y ) i[8,t] !i )u;{a,t) x ^. 
The second term is treated exactly in the same way. Summarizing, we have shown that 

CyV 7 (i? (x _ y) ;M 2 MM)^- 



t p-t 

X 8 , u (g> dX u - / y. u ® dY u 



< 



L 2 



Finally, by Lemma 13.21 

and by the Cauchy-Schwarz inequality 



Voo(R{x-Y);[s,t] 2 ) <4 sup \X U -Y U \\ 2 
ue[s,t] 



which gives the claim. 



□ 



Corollary 3.7. Under the assumptions of Provosition W7b\ for every 7 satisfying 7 > p and 
- + ->!, and every p > 2p and 7' > 7, there is a (random) control ui such that 



|X» t |<£( S ,t) n/p 
|Y s " )t |<£( S ,*) n/p 



\Tq it -Y? tt \<e&( 8 ,t)W+— 

holds a.s. for all s < t and n — 1,2 where e = sup uG [ T i \X U 
is a constant C — C(p, p, j,y,K) such that 



(4) 
(5) 
(6) 

Y u \jji"' ■ Furthermore, there 



\cj(0,T)\ Lg <CT(qP/ 2 + q^) 

holds for all q > 1 . 

Proof. Let qj be the control from Proposition 13. 61 We know that 

\Xl t -Yl t \ L2 < Cl eoj{ Sl t)^+^ 

holds for a constant c\ for all s < t and n = 1,2. Furthermore, |X™ t | L 2 < C2w(s,t)2p for a 
constant C2 for all s < t and n = 1, 2 and the same holds for Y (this just follows from setting 
Y = const, and 7 = p in Proposition ^. 61) . Now introduce a new process X : [0, T] — >• R d on 
the same sample space as X such that for all sample points, we have 

X w (o,t)/ w (o,T) = Xt, Vt e [0,T], 

and define Y in the same way. Then X, Y are well defined, multiplicative, and we can 
replace the control ui by c^K\t — s\ for the two re-parametrized processes. Using that X, Y 



The reader might feel a bit uncomfortable at this point asking why it is allowed to put expectation inside 
the integral (which is not even an integral in Riemann-Stieltjes sense). However, this can be made rigorous by 
dealing with processes which have sample paths of bounded variation first and passing to the limit afterwards 
(cf. [FVlOal iFVlObl IFR121 IFH12] 1. We decided not to go too much into detail here in order not to distract the 
reader from the main ideas and to improve the readability. 
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= n™/2 



are Gaussian, we may pass from L to L q estimates and we know that 0(|X™|l<j) = q 
(same for Y and X - Y, cf. [FVlObl Appendix A]). Hence 



\Xl t \ L «<c4VqKVp) n \t-s\% 



\Yl t \ Lq <c 4 (VqIO/Pr\t~s\^ 



ft - ft| £ , < I Ci (VqKVPr\t - s\^+^r 



(7) 
(8) 
(9) 



hold for all s < t, n = 1,2 and q > 1 with e = eif 5 ^ 2 p . Using Lemma 16. II in the appendix, 
we see that there is a constant C5 = Cs(p, p, 7, 7', if) such that 



ft I 



sup 

S <*G[0,T] I* " 



SUp : 7—r 

«te[o,T] \t - s\"/p 



L" 



Li 



sup 



ft-Yy t | 



L9 



s <*e[o,T] It-spMn) 
hold for g sufficiently large and n — 1,2 where = ^7 + ^^T"- Set 

™p J2\*u,nJ p/n 

DC[.,t] tiGjD 

«?(*,t):= sup El Y ".,* I+1 l P/ " 

DC[«,t] tiGjD 

£_ y (M):= sup £ |X£ )ti+i -Y™ i4+ J^ 

DC[«,t] { . eD 



(10) 

(11) 
(12) 



and 



£>(«,<):= Wx(s,i)+w?(s,t) + e5Wwl_ y (s,i). 



n=l,2 



for s < t. Clearly, w fulfils (Q}, (|S|) and ©. Moreover, the notion of p- variation is invariant 
under reparametrization, hence 

lx™ \p/ n 

*£(0,T)= sup £ |X» itt+1 |f/"= sup ^ft ti+i K«<T sup i-^— 

and a similar estimate holds for d)y(0,T) and a)^_ y (0,T). By the triangle inequality and 
the estimates ©, © and ©, 

|^(0,T)| Lg < 2 |o>^(0,T)| i , + |a>?(0,T)| i ,+epfep5_ y (0,T)| z , 

n=l,2 

< c 6 T (qP/ 2 +q E ¥+ q^ 2 )) < C 7 T(q*>/ 2 + q<) 

for q large enough. We can extend the estimate to all q > 1 by making the constant larger 
if necessary. □ 



Corollary 3.8. Let Co fee the random control defined in the previous corollary. Then, for 
every n, there exists a constant c n such that 

|X" t | < c n Co(s,t)f, |Y? )t | < Cj&fat)? 

a.s. for all s < t. The constants c n are deterministic and can fee chosen such that c n < ^/ p y. > 
where x\ :— T(x — 1). 

Proof. Follows from the extension theorem, cf. |L98[ Theorem 2.2.1] or [LCL06J Theorem 
3.71. □ 
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4 Main estimates 

In what follows, we let p E (2p, 3). Let 7 > p such that - + - > 1. We write log + x — 
max{x,0}, and set 

e= sup |X u -r u | L2 T 

u£[0,T] 



4.1 Higher level estimates 

We first introduce some notations. Suppose X is a multiplicative functional in T N (R d ) 
with finite p- variation controlled by ui, N > \p\. Then, define 

N 

X S)i = 1 + X M E T N+1 (R d ). 

n=l 

Then, X is multiplicative in T , but in general not in T N+1 . For any partition D = {s = 
u < ui < ■ ■ ■ < ul < ul+i = t}, define 

X£ t := X SjUl <g> •• • <g> X„ i)t G T^+^M 4 ). 

The following lemma gives a construction of the unique multiplicative extension of X to 
higher degrees. It was first proved in Theorem 2.2.1 in L98 . 

Lemma 4.1. Let X be a multiplicative functional in T . Let D = {s < u\ < ■ ■ ■ < mj, < t] 

be any partition of (s,t), and D J denote the partition with the point Uj removed from D. 
Then, 

N 

X£ t - Xg = K-r* ® G T N +\R d ). (13) 

n=l 

In particular, its projection onto the subspace T N is the 0-vector. Suppose further that X 
has finite p-variation controlled by uj, and N > [pj , then the limit 

lim ±f t e T N+1 (R d ) 

exists. Furthermore, it is the unique multiplicative extension of X to T N+1 with finite p- 
variation controlled by ui. 

Theorem 4.2. Let (X, Y) and p, 7 as in Provosition lS. 61 Then for every p > 2p and 7' > 7 
there exists a constant C3 depending on p and 7' and a (random) control uj such that for all 
q>l, we have 

|(&(0,r)|i« < M < +00, 
where M — M{p, ^,7,7' , K,q), and the following holds a.s. for all [s,t]: 
(i) V W + § > I, then 

\Xl t -Yl t \<C 3 euj(s,t)W+ 1 P . 

(11) If ^ + f = l, then 

|X* t - Y% t \ < C 3 e ■ (1 + log+ [«(Q, T)/e 1 -^] ) • oj(s, t). 
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(Hi) If 2y + f < 1, i/ien 

|X 3 t -Y 3 { | <C 3 eray^( s ,0, 

Proof. Let s < t E [0, T] and let cD be the (random) control defined in Corollary 13.71 Then, 
by the same corollary, for every q > 1, \oj(0, T)\li < M. Fix an enhanced sample rough path 
(X, Y) up to level 2 and for simplicity, we will use uj to denote the corresponding realisation 
of the (random) control uj. We can assume without loss of generality that 

e < u(s,*)*~^, (14) 

otherwise there will be nothing to prove. Let D = {s = Uq < ■ ■ ■ < ul + \ = t} be a 
dissection. Then (cf. |L98[ Lemma 2.2.1]), there exists a j such that 

2 

uj(uj-i,Uj + i) < — u)(s,t), L > 1. (15) 
Let -D- 7 denote the dissection with the point Uj removed from D. Then, we have 

|(X£ t - Y£) 3 | < |(xg - Y s ^) 3 | + E(I<_ 1 ,„, ®x^ 3 . +1 | 

k=l 

where R s t = Y Sjt — X Sj t. By assumption, 

® R'7X +1 1 < C • min |e(|o;( S , t)) ^ , (|w( S , t)) 'V (16) 
and similar inequalities hold for the other two terms in the bracket. Thus, we have 

|(Xf t -Y s c t ) 3 | < -Y^) 3 | + C 3 min| e (i W ( S ,t))^^,(i W ( S ,t))t|. 

Let iV be the integer that 

{-L-ujis,^-^ < e < [±u(s,tp-W, (17) 

then 

e[-uj(s,t)]^' + p < [~w(s,i)]» 

if and only if L < iV. By Lemma T4. 11 we have 

X 3 t = lim (Xf t ) 3 , Y s 3 t = lim (Yf t ) 3 . 

Thus, for a fixed partition D, we choose a point each time according to (|15[) . and drop them 
successively. By letting \D\ — > +oo, we have 



X 3 . t - Y 3 1 | < C 3 



AT 1 +oo 



L=l L=7V+1 



Approximating the sums by integrals, we have 

- N , , ,, „ /•+<» 



|Xf, t 



1 2 /* 1 2 f^r°° 

Y s.t\ < C 3 [eoj(s,t)^ T + p(l + / x _( ^ 7 ' + p ) di) + a;(s,*)t / aT*da;]. 

Ji Jn 

Compute the second integral, and use 



l N+ 1 V ' n 



we obtain 



|X 3 t -Y 3 t | < C 3 [euj(s,t)^ T + i(l + j x~ ( ^ T+ ^ ) dx) + e T ^ r uj(s,t)}. (18) 
Now we apply the above estimates to the three situations respectively. 
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1. A + - > 1. 

27 p 
In this case, the integral 



N r+00 



'1 Ji 
converges. On the other hand, (fTl)) implies 

3-p 1 I 2 

ei-p/Vo)(s,t) < ew(s,i) 5 T 7 

thus, from (fT51) . we get 

|X* -Y» |<C 3 e W ( s ,t)^ + i 



2. -L. + 2 = 1. 

In this case, | — ^7 = and = 1- Thus, by the second inequality in (|17[) . wc 

have 

/ x dx — log TV < log w(s,t) loge. 

J i 3-p 

On the other hand, (fT^)) gives 

log w(s,t) - — — loge > 0. 
3-p 

Combining the previous two bounds with (|18l) . we get 

|XL - Y s 3 t | < C 3 e[l+loga;( Sj t) _ _£_]og e ] w ( a ,t). 

3-p 

We can simplify the above inequality to 

|X* t - Y%\ < C 3 e[l + log + (u(0,T)/e^)Ms,t), 

where we have also included the possibility of e > w(0,T)p _1 . 
Now we have 

I 1 1 I 2 \ 1 1 2 _r-t 1 2wl 1_ 1 1 2 

1 + / x ( *r'^p'dx < CN 27' p<C-e { 27' Vwfs,;) 1 27' p, 



where the second inequality follows from (1171) . Combining the above bound with (|18p . wc 
obtain 



\Xl t -Yl t \<C 3 e^7Wu(s,t). 



□ 



The following theorem, obtained with the standard induction argument, gives estimates 
for all levels n = 1, 2, • • • . 
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Theorem 4.3. Let (X,Y) and p, 7 as in Provosition \3. 61 p > 2p and 7' > 7. TTien i/iere 
exists a (random) control u) such that for every q > 1, we have 

P(0,T)| ia < M 

where M = M (p, p, 7, 7', g, K), and for each n there exists a (deterministic) constant C n 
depending on p and 7' such that a.s. for all [s,t\: 

(i) If 2y + I > 1} i/ien we /iai>e 

|X? )t -Y? )t | <C n e<2;(M)^ + ^ 
(lij J/ + - = 1> then we have 

|X« t - Y« t | < C„e- (l + log+ [ifO.TJ/e 1 -*]) -a;( S ,t)^ + ^. 
fm,) J/ + I < 1, i/ien /or all s < t and all small e, we have 

|X:. t -Y^<C„ e T^^ (M) ^™. (19) 

Proof. We prove the case when ^7 + - < 1; the other two situations are similar. Let u> be 
the control in the previous theorem. Fix an enhanced sample path (X, Y), the corresponding 
realisation ui of &, and s < t G [0, T]. We may still assume (|T4"f without loss of generality. 
Thus, for n = 1, 2, we have 

|X" t - Y£ t | < eo;(s,t)^ 7+:i ? L < C^ra^s, i) " ^ {P \ 

where the second inequality comes from (I14[) . The above inequality also holds for fc = 3 by 
the previous theorem. Now, suppose (|19[) holds for k — 1, • ■ ■ , n, where n > 3, then for level 
fc = n + 1, the exponent is expected to be 

n + {p} 

so that the usual induction procedure works (cf. |L98] . Theorem 2.2.2.). Thus, we prove 
(3D for all n. □ 



4.2 Proof of Theorem [TTT1 

Proof. We prove the second situation when ^- + - < 1. The first one is similar. Let 

e = sup^gjo t] l-^u — ^"1^2 ~' ■ It is sufficient to show that for every p > 2p there is a constant 
C such that 

JV 3 -UL- 



|^_ var (X,Y)| i5 <Ce 



1— p/i 



where a > 2-f and N > [a\ both satisfy the assumptions of Theorem 11.11 Set 

p' :=(l+r))p, p:=2{l + 2 V )p, 7' := (1 + 77)7, 7" := (1 + 277)7 

for some 77 > 0. We can choose 77 small enough such that y + y > 1 and p < 3 hold, and 
the conditions of Theorem 14.31 are satisfied for p' and 7'. Clearly yr + -<^- + ^<l, thus 

|X" - Y" I < C n e^hu{s,t)'^ iSl 



Theorem 14.31 implies that 



holds a.s. for any n and s < t where ui is a random control as in Theorem l4.3l Furthermore, 
for any n, 

p 2 7 " p 2 7 " ; p 2 7 " v n p 2i'> v p 27" 7 

(20) 
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Note that the last expression implies that 

^:=n(i-JL)-^M>0 
p 2y' p 

for all n. Fix a dissection D = {0 = uo < . . . < ul < T} of the interval [0, T]. Using 
ui(ui,Ui+i) < w(0, T), we have 

z ' ^ i 

Choosing 77 smaller if necessary, we may assume that a > 27" and super-additivity of the 
control implies 



-i+{p} 



■ i * 

Passing to the supremum over all partitions of [0,T], we have 

sup (j2 \*l,u i+1 ~ Vl, Ui+1 \A ' < C n e^h^,Tf 

D \ i / 

Let q > 1. By Thcorcm l4.3[ there is a constant M depending on p, 7, a, S, q and K such that 
|cD(0,T)|i« < M. Taking L 9 norm on both sides, we have 

|^_ var (X,Y)| L9 <Ce^ 
which was the claim. □ 



5 Applications 

5.1 Convergence rates of rough differential equation 

Consider the rough differential equation of the form 

d 

dY t = Vi{Y t ) dX\ =: V{Y t ) dX t ; Y e M. e (21) 
i=i 

where X is a centered Gaussian process in ~R d with independent components and V = {Vi)f =1 
a collection of bounded, smooth vector fields with bounded derivatives in K. e . Rough path 
theory gives meaning to the pathwise solution to (|2T|) in the case when the covariance Rx 
has finite p- variation for some p < 2. Assume that p £ [1, | ) and that there is a constant K 
such that 

y p (i?x;M 2 ) < K\t-s\ k p (22) 

for all s < t (note that this condition implies that the sample paths of X are a-H61der for 
all a < j^). For simplicity, we also assume that [0, T] = [0, 1]. For every k £ N, we can 
approximate the sample paths of X piecewise linear at the time points {0 < 1/k < 2/k < 
. . . < (k — l)/k < 1}. We will denote this process by X^ k \ Clearly, X^ — !> X uniformly 
as k — > 00. Now we substitute X by X^ in (|2"Tj) . solve the equation and obtain a solution 
Y W ; we call this the Wong-Zakai approximation of Y. One can show, using rough path 
theory, that Y^ -4 Y a.s. in uniform topology as k — > 00. The proposition below is an 
immediate consequence of Theorem 11.11 and gives us rates of convergence. 
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Proposition 5.1. The mesh size 4 Wong-Zakai approximation converges uniformly to the 
solution of (|21jl with a.s. rate at least fc^ 5 p~ 1 ~ <5 - ) for any 6 S (0, — 1). In particular, the 
rate is arbitrarily close to ^ urfien p = 1, which is the sharp rate in that case. 

Proof. First, one shows that (|22p implies that 

sup \x[ k) -X t \ L 2=0{k-^). 
te[o,i] 

One can show (cf. jFVlObl Chapter 15.2.3]) that there is a constant C such that 

supFpORp^w^M 2 ) < C\t-s\p 
fceN 

holds for all s < t. By choosing g large enough, a Borel-Cantelli type argument applied to 
Theorem 11.11 shows that g CT _ var (X, X' fc ^) — > a.s. for k — > oc with rate arbitrarily close to 

, -iri-i-i .„ 1 1 

fc 21 p ->>, if h - > 1, 

2 7 p 

and arbitrarily close to 

k-^\ if — + - < 1, 
2 7 p 

both cases are subject to 7 > | and - + - > 1. Note that in the second situation, the 

actual value of 7 does not matter, and we always have a rate of 'almost' j — 1. For the 
first situation, we need to let 7 as large as possible but still satisfy the constraints. The 
critical value is = ^(1 — -), which also results in a rate that is arbitrarily close to j- — 1. 
Using the local Lipschitz property of the Ito Lyons map (cf. |FV10b[ Theorem 10.26]), we 
conclude that the Wong-Zakai convergence rate is faster than 

for any 8 > (but not for <5 = 0). □ 

Remark 5.2. For p G (1, §), the rate above is not optimal. In fact, the sharp rate in this 
case is 'almost' - — - \, as shown in WRIS^ . The reason for the non-optimality of the rate 
is that we obtain the third level estimate merely based on the first two levels, which leads to 
a reduction in the exponent in the rate. On the other hand, this method does not use any 
Gaussian structure on the third level, and can be applied to more general processes. For the 
case p = I, we recover the sharp rate of 'almost' 5. 



5.2 The stochastic heat equation 

In the theory of stochastic partial differential equations (SPDEs), one typically considers 
the SPDE as an evolution equation in a function space. When it comes to the question 
of time and space regularity of the solution, one discovers that they will depend on the 
particular choice of this space. As a rule of thumb, the smaller the space, the lower the time 
regularity •: III)!) . Section 5.1). The most prominent examples of such spaces are Hilbert 
spaces, typically Sobolev spaces. However, in some cases, it can be useful to choose rough 
paths spaces instead f [Hll] ]). A natural question now is whether the known regularity results 
for Hilbert spaces are also true for rough paths spaces. In this section, we study the example 
of a modified stochastic heat equation for which we can give a positive answer. 

Consider the stochastic heat equation: 

dtp = (d xx - l)ip dt + adW (23) 
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where a is a positive constant, the spatial variable x takes values in [0, 2n], W is space-time 
white noise, i.e. a standard cylindrical Wiener process on L 2 ([0, 27r],]R d ), and ip denotes 
the stationary solution with values in M. d . The solution t/j is expected to be almost 4- 
Holder continuous in time and almost ^-Holder continuous in space (cf. H09 ). In the next 
Theorem, we show that this is indeed the case if we choose the appropriate rough paths 
space. 



Theorem 5.3. Let p > 2. Then, for any fixed t > 0, the process x <— > ipt(x) is a Gaus- 
sian process (in space) which can be lifted to an enhanced Gaussian process ^t(')> a pro- 
cess with sample paths m C°' p -™ r ([0, 2tt], G^ (R d )). Moreover, t i-> * t (-) has a Bolder 
continuous modification (which we denote by the same symbol). More precisely, for every 

a G K), \ — 2~p)> th ere exists a (random) constant C such that 

Qp- va r(y s ,Vt) <C\t- S \ a 

holds almost surely for all s < t. In particular, choosing p large gives a time regularity of 
almost j-Holder. 

Proof. The fact that x i— ► ipt ( x ) can be lifted to a process with rough sample paths and that 
there is some Holder-continuity in time was shown in Lemma 3.1 in [Hllj . see also [FH12 . 
We quickly repeat the argument and show where we can use our results in order to derive 
the exact Holder exponents. Using the standard Fourier basis 



•jm sin(fciz;) if k > 







i- if k = 

Z7T 



e k (x) = 

^ cos(fcx) if k < 

the equation (|23D can be rewritten as a system of SDEs 

dY t k = -{k 2 + l)Y t k dt + a dW t k 

where (W k )kez is a collection of independent standard Brownian motions and (Y k )kez are 
the stationary solutions of the SDEs, i.e. a collection of centered, independent, stationary 
Ornstein-Uhlenbeck processes. The solution of (|23[) is thus given by the infinite sum ipt{%) — 
J2kez Yt ke k( x )- One can easily see that 

where Id denotes the identity matrix in R dxd . In particular, for s = t, 

E [ip t (x) (8) My)} = K(x - y) x I d 

where K is given by 

K(x) = , . , . . cosh(|a;| — n) 
4smh(7r) 

for x G [— 7r, 7r] and extended periodically for the remaining values of x (this can be derived 
by a Fourier expansion of the function x i— > cosh(|x| — n)). In particular, one can calculate 
that x i — y i^tix) is a Gaussian process with covariance of finite 1-variation (see the remark at 
the end of the section for this fact), hence ipt can be lifted to process ^>t with sample paths 
in the rough paths space C°< p - var ([0, 2tt], G*LpJ (R d )) for any p > 2. 

Furthermore, for any s < t, x t-> {ip s {x),^)t(x)) is a Gaussian process which fulfils the 
assumptions of Theorem ll.ll and the covariance Ri^^^a also has finite 1-variation, uniformly 
bounded for all s <t, hence 

SUp |-R(^ s ,^ t )|l-var;[0,27r] 2 = : c l < °0- 
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Therefore, for any 7 £ (l,p/2) and g > 1 there is a constant C = C(p, 7, ci, g) such that 

\Qp- vaI (^ s , *t)li9 < C SU P FtW - V' S (^)I L 2 7 
:re[0,27r] 

holds for all s < t. A straightforward calculation (cf. |H111 Lemma 3.1]) shows that 

\Mx)-Mx)\l» <c 2 \t-s\ 1/4 
for a constant ci. In particular, we can find 7 and g large enough such that 

Ki-^-i = A_J_A_i < i_J_ 

9 \4 47/ g 4 2p' 

Since C°' p ~ vaI is a Polish space, we can apply the usual Kolmogorov continuity criterion to 
conclude. □ 



Remark 5.4. We emphasize that here, for every fixed t, the process ipt(-) is a Gaussian 
process, where the spatial variable x should now be viewed as 'time'. This idea is due to 
M.Hairer. Knowing that the spatial regularity is 'almost' 1/2 for every fixed time t, one 
could guess that covariance of this spatial Gaussian process has finite 1-variation. For a 
formal calculation, we refer to \H09jj or \FH1^ . 



6 Appendix 



The next Lemma is a slight modification of |FV10b[ Theorem A. 13]. The proof follows the 
ideas of [FH121 Theorem 3.1]. 

Lemma 6.1 (Kolmogorov for multiplicative functional) . Let X,Y: [0,T] x Q -> T N (V) 
be random multiplicative functionals and assume that X(w) and Y(w) are continuous for all 
u> £ CI. Let j3, 5 £ (0,1] and choose /3' < j3 and 5' < 8. Assume that there is a constant 
M > such that 



s,t\Ll/' 



<M n \t-s 
< M n \t 



nd 



m/3 



|X? t -Y? t | is/ n<M"e|t- a 



hold for all s < t £ [0, T] and n = 1, . . . , TV where e is a positive constant and q > qo where 

1 1 



9o 



1 



V 



J-P' 8-8' 
Then there is a constant C = C{N,j3,j3' ,8,8') such that 



|X" t | 



s<te 



sup 

s<te[o,T] 



\yi\t\ 



\nj3' 



|X?, t -Y? |f | 

S <*6[0,T] \t-s\S' + (n-DP> 



< CM n 



< CM n 



< CM n e 



(24) 
(25) 
(26) 



hold for all n = 1, . . . , N. 
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Proof. W.l.o.g., we may assume T = 1. Let (Dk)ken be the sequence of dyadic partitions 
of the interval [0, 1), i.e. D k = { ^ : I = 0, . . . , 2 k - l}. Clearly, \D k \ = ^ = 2~ k . Set 



Kl x := max X™. t 
Kl Y := max |Y t n t \ 

K lx-Y :=\™g.\XZ, u+1 -YZ tti+1 \ 



for n = 1, . . . , N and fc 6 N. By assumption, we have 

E\Kl x \i <EJ2 l X *U +1 l" < #D k max B|X^ ti+ J* < M^^" 1 . 

ti£D k 

In the same way one estimates K^ Y an d K k x-Yi nence 

\K£ x \ Lq/n < M n \D k \ nf> - n/q (27) 
\Kl Y \ Lq/n < M n \D k \ n P- n /« (28) 
|Xfc,x-yL«/» ^M"!^!^- 1 ^-"/". (29) 

Note the following fact: For any dyadic rationals s < t, i.e. s < t £ A := [J^L-l -Dfc, there is 
araeN such that |Z? m +i| < \t — s\ < \D m \ and a partition 

s = t < ti < . . . < tat = t (30) 

of the interval [s,t) with the property that for any i = 0, . . . , N — 1 there is a k > m + 1 
with [ri, Tj+i) G -Dfc, but for fixed fc > m + 1 there are at most two such intervals contained 
in D k . 

Step 1: We claim that for every n = 1, . . . , N there is a real random variable K x such 
that \K x \ Lg / n < M n c where c = c{(3, (3' ,5,5') and that for any dyadic rationals s < t and 



m, (Ti)^L chosen as in (|3"0l we have 



N-l l-v-ra I 

y ^ < (3i) 

z=0 1 1 

Furthermore, the estimate (|31l) also holds for Y" and a random variable K Y . Indeed: By 
the choice of m and {n)fL , 

N—l \~v-n I oo 0T^ n 00 J^ n 00 

E I A r,-,T i + il \ - ZIX k,X <9 \ - fVfc^f „ „ ra 

| t _ s |«/3< - 2-, lAn+lK - \ 2 1 1 IW' " fl l^l""' X ' 
2=0 fc=m+l k=m+l fc=l 

It remains to prove that -ff£| L ,,/„ < M"c. By the triangle inequality and the estimate (|2"T|) . 



E 



fc=i 1 1 



<M"^ |D fe |"( /3 - 1 /«-' 3 ') < M" y \D k \^- 1 l q °- f3 '^ < oo 

fe=l fc=l 



since /3 — 1/qo — ft' > which shows the claim. 



Step 2: We show that (1241) and ([25]) hold for all n = 1, . . . , N. It is enough to consider X. 

IX 11 I 

Note first that, due to continuity, it is enough to show the estimate for sup s<tgA . . 
By induction over n: For n = 1, this just follows from the usual Kolmogorov continuity 
criterion. Assume that the estimate is proven up to level n — 1. Let s < t be any dyadic 
rationals and choose m and (t^^Lq as in (|30[) . Since X is a multiplicative functional, 

N-l n-l N-l 

ixj,«i < 2 i^+J+E^J*^ 1 E i x U+J 

i=0 1=1 ' ' i=Q 
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and thus, using step 1, 



Xn I N—l lyn I n-1 |vn-l| AT-1 i-y-i I 

s,t\ < \ "* l A Tj.T i + i \ - ^ I \ - l A n,r i+ il 

| t _ s |n/3< " | t - a |nfl' + Z^ u ^P A | 1; _ u |(n-0^ 2v | t - «|W 

n—l lY«-'l 



<^+x; bu P ' 



( 



We can now take the supremum over all s < t S A on the left. Taking the L 9 /™-norm on 
both sides, using first the triangle, then the Holder inequality and the estimates from step 
1 together with the induction hypothesis gives the claim. 

Step 3: As in step 1, we claim that for any n = 1, . . . , N there is a random variable 
K%_ Y e Lil n such that for any dyadic rationals s < t and to, (Ti)fL chosen as above we 
have 

N—l ivji -yn 

E < 02) 

i— ' ' 

Furthermore, we claim that I-Kjc-yI-L"/™ < M n c where c = c(/3, j3', 8, 5'). The proof follows 
the lines of step 1, setting 

N—l Ivn -\rn I oo j^n 



f \l _ s |<5'+(n-l)/3' - Z Z-j |£) fe |5'+(ra-l)/9' ' - 



Step 4: We prove that (J2S]) holds for all n = 1, . . . , N. By induction over n: The case 
n = 1 is again just the usual Kolmogorov continuity criterion applied to t h- > e (JQ — if). 
Assume the assertion is shown up to level n—l and chose two dyadic rationals s < t. Using 
the multiplicative property, we have 



N-i n-1 N-i 



\Kt ~ Y "*l ^ E l X ",r, +1 - Y^, Ti+1 | + Z ^ i max^lX^I £ |X^ n+1 - Y^ Ti+1 1 

i=0 (=1 i=0 

n-1 N-1 

+ E i= T* N l X "n - Y "n I E l Y t 4 ,r i+1 l- 



Now we proceed as in step 2, using the estimates from step 1 to step 3 and the induction 
hypothesis. □ 
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